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Abstract. The spacetime geometry on the equatorial slice through a Kerr black 
hole is formally equivalent to the geometry felt by phonons entrained in a rotating 
fluid vortex. We analyse this example of "analogue gravity" in some detail: First, we 
find the most general acoustic geometry compatible with the fluid dynamic equations 
in a collapsing/expanding perfect-fluid line vortex. Second, we demonstrate that there 
is a suitable choice of coordinates on the equatorial slice through a Kerr black hole 
that puts it into this vortex form; though it is not possible to put the entire Kerr 
spacetime into perfect-fluid "acoustic" form. Finally, we discuss the implications of 
this formal equivalence; both with respect to gaining insight into the Kerr spacetime 
and with respect to possible vortex- inspired experiments, and indicate ways in which 
more general analogue spacetimes might be constructed. 



PACS numbers: 02.40.Ma. 04.20.Cv, 04.20.Gz, 04.70.-s 
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1. Introduction 



It is by now well-known that the propagation of sound in a moving fluid can be described 
in terms of an effective spacetime geometry. (See, for instance, [1, 2, 3, 4, 5, 6], and 
references therein). 

• In the geometrical acoustics approximation, this emerges in a straightforward 
manner by considering the way the "sound cones" are dragged along by the fluid 
flow, thereby obtaining the conformal class of metrics (see, for instance. 



-Vi 



-Vi 



hi 



(1) 



Here c is the velocity of sound, v is the velocity of the fluid, and hij is the 3-metric 
of the ordinary Euclidean space of Newtonian mechanics (possibly in curvilinear 
coordinates). 

In the physical acoustics approximation we can go somewhat further: A wave- 
equation for sound can be derived by linearizing and combining the Euler equation, 
the continuity equation, and a barotropic equation of state [1, 2, 3, 4]. For 
irrotational flow this process leads to the curved-space d'Alembertian equation and 
in particular now flxes the overall conformal factor. The resulting acoustic metric 
in (3-1-1) dimensions is (sec, for instance, [4, 5]): 
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Here p is the density of the fluid. Sound is then described by a massless minimally 
coupled scalar fleld propagating in this acoustic geometry. In the presence of 
vorticity a more comphcated wave equation may still be derived [7] . For frequencies 
large compared to the local vorticity u = ||V x ■i7|| that more complicated wave 
equation reduces to the d'Alembertian and the acoustic geometry can be identifled 
in the usual manner. For more details see [7]. 

Now because the ordinary Euclidean space {hij) appearing in these perfect fluid acoustic 
geometries is Riemann flat, the 3-dimensional space (given by the constant-time slices) in 
any acoustic geometry is forced to be conformally flat, with 3-mctric Qij = (p/c) hij. This 
constraint places very strong restrictions on the class of (3+l)-dimensional geometries 
that can be cast into perfect fluid acoustic form. While many of the spacetime 
geometries of interest in general relativity can be put into this acoustic form, there 
are also important geometries that cannot be put into this form; at least for a literal 
interpretation in terms of flowing perfect fluid liquids and physical sound. 

In particular, the Schwarzschild geometry can, at the level of geometrical acoustics, 
certainly be cast into this perfect fluid acoustic form [4]. However, at the level of physical 
acoustics (and working within the context of the Painleve-Gullstrand coordinates) 
there is a technical difficulty in that the Euler and continuity equations applied to 
the background fluid flow yield a nontrivial and unwanted (3-1-1) conformal factor 
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p/ c [4] . This overall conformal factor is annoying, but it does not affect the sound cones, 
and so does not affect the "causal structure" of the resulting analogue spacetime [6]. 
Furthermore any overall conformal factor will not affect the surface gravity [8], so it will 
not directly affect the acoustic analogue of Hawking radiation [1,4]. So for most of the 
interesting questions one might ask, the Schwarzschild geometry can for all practical 
purposes be cast into acoustic form. 

Alternatively, as we shall demonstrate below, the Schwarzschild geometry can also 
be put into acoustic form in a different manner by using isotropic coordinates — see the 
appendix of this article for details. More generally, any spherically symmetric geometry 
(static or otherwise) can be put into acoustic form — details are again provided in the 
appendix. (This of course implies that the Reissner-Nordstrom geometry can be put 
into acoustic form.) In particular, the FRW cosmologies can also be put into acoustic 
form, both at the level of physical acoustics and at the level of geometrical acoustics. 
In fact there are two rather different routes for doing so: Either by causing the fluid to 
explode or by adjusting the speed of sound [9, 10, 11, 12, 13, 14, 15, 16]. 

There is however a fundamental restriction preventing the Kerr geometry [and 
Kerr-Newman geometry] being cast into perfect fluid acoustic form. It has recently 
been established [17, 18, 19] that no possible time-slicing of the full Kerr geometry 
can ever lead to conformally flat spatial 3-slices. Faced with this fact, we ask a more 
modest question: Can we at least cast a subspace of the Kerr geometry into perfect fluid 
acoustic form? Specifically, since we know that the effective geometry of a generalized 
line vortex (a "draining bathtub" geometry) contains both horizons and ergosurfaces 
[4, 5], one is prompted to ask: If we look at the equatorial slice of the Kerr spacetime 
can we at least put that into acoustic form? If so, then this opens the possibility of 
finding a physically reasonable analogue model based on a vortex geometry that might 
mimic this important aspect of the Kerr geometry. Thus we have three independent 
physics questions to answer: 

• What is the most general perfect fiuid acoustic metric that can [even in principle] 
be constructed for the most general [translation invariant] line vortex geometry? 

• Can the equatorial slice of Kerr then be put into this form? 
(And if not, how close can one get?) 

• By generalizing the analogue model to something more complicated than a perfect 
fiuid, can we do any better? 

We shall now explore these three issues in some detail. 

2. Vortex flow 

2. 1 . General framework 

The background fiuid fiow [on which the sound waves are imposed] is governed by three 
key equations: The continuity equation, the Euler equation, and a barotropic equation 
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of state: 



|^ + V-(p^;) = 0. (3) 



- + (..V). 



-Vp + /. (4) 
p = p{p). (5) 



Here we have included for generality an arbitrary external force /, possibly magneto- 
hydrodynamic in origin, that we can in principle think of imposing on the fluid flow to 
shape it in some desired fashion. From an engineering perspective the Euler equation is 
best rearranged as 

-Vvp, (6) 



/ 



with the physical interpretation being that / is now telling you what external force you 
would need in order to set up a specified fluid flow. 

2.2. Zero radial flow 

Assuming now a cylindrically symmetric time-independent fluid flow without any sinks 
or sources we have a line vortex aligned along the z axis with fluid velocity v: 

v{r) = Vg{r) 9. (7) 

The continuity equation (3) for this geometry is trivially satisfied and calculating the 
fluid acceleration leads to 

U^(v.V)v^-^^^f. (8) 

Substituting this into the rearranged Euler equation (6) gives 

/ = ff r = |-p(r) + c\r) a,p(r)| f (9) 

with the physical interpretation that: 

• The external force ff must be chosen to precisely cancel against the combined 
effects of centripetal acceleration and pressure gradient. The angular- flow is 
not completely controlled by this external force, but is instead an independently 
specifiable quantity. (There is only one relationship between /^(r), p(r), c(r), and 
ve{r), which leaves three of these quantities as arbitrarily specifiable functions.) 

• We are now considering the equation of state to be an output from the problem, 
rather than an input to the problem. If for instance p(r) and c^(r) are specified 
then the pressure can be evaluated from 

p{r) = j c\r) ^ dr, (10) 

and then by eliminating r between p{r) and p(r), the EOS p{p) can in principle be 
determined. 
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For zero external force (and no radial flow), which arguably is the most natural 
system to set up in the laboratory, we have 

vl{r) 



r dr In p(r), (11) 



which still has two arbitrarily specifiable functions. 

In the geometric acoustics regime the acoustic line-element for this zero-source/ 
zero-sink line vortex is 

ds^ oc - (c^ - vl) dt^ - 2vg rdedt + dr^ + r^e^ + dz\ (12) 

In the physical acoustics regime the acoustic line-element is 

ds^ = - (c^ - vj) dt^ - 2vq rde dt + dr^ + r'^dO^ + dz"^ . (13) 

The vortex quite naturally has a ergosurface where the speed of the fluid flow 
equals the speed of sound in the fluid. This vortex geometry may or may not have 
a horizon — since v.f is identically zero the occurrence or otherwise of a horizon 
depends on whether or not the speed of sound c(r) exhibits a zero. X We shall later 
see that this class of acoustic geometries is the most natural for building analogue 
models of the equatorial slice of the Kerr geometry. 



2.3. General analysis with radial flow 

For completeness we now consider the situation where the vortex contains a sink or 
source at the origin. (A concrete example might be the "draining bathtub" geometry 
where fluid is systematically extracted from a drain located at the centre.) Assuming 
now a cylindrically symmetric time-independent fluid flow with a line vortex aligned 
along the z axis, the fluid velocity v is 

V = Vf.{r) r + Vfj{r) 9. (14) 

Wherever the radial velocity Vr{r) is nonzero the entire vortex should be thought of as 
collapsing or expanding. 

The continuity equation (3) for this cylindrically symmetric problem is 

V-(pt?) = 0, (15) 

and the rearranged Euler equation (6) for a pressure p which depends only on the 
radial-coordinate r is 

f ^ p(v-V)v + c^ drpf. (16) 

I The vanishing of c(r) at a horizon is exactly what happens for Schwarzschild black holes (or their 
analogues) in cither Schwarzschild or isotropic coordinate systems. This is rather different from the 
behaviour in Painleve-GuUstrand coordinates, but is a quite standard signal for the presence of a 
horizon. 
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Figure 1. A simple non-collapsing vortex geometry: The green circles denote the 
streamlines of the fluid Row. The outer red circle represents the ergosurface, where the 
fluid velocity reaches the speed of sound. The inner red circle (which need not exist in 
general) represents a possible [outer] horizon where the speed of sound goes to zero. 




Figure 2. A collapsing vortex geometry (draining bathtub): The green spirals denote 
streamlines of the fluid flow. The outer circle represents the ergosurface while the inner 
circle represents the [outer] event horizon. 

At this stage we note that p is in general not an independent variable. Because equation 
(15) corresponds to a divergence-free field, integration over any closed circle in the two- 
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dimensional plane yields 

j) p{r) v{r) ■ f ds = 27r p{r) Vf{r) r = 27r ki. 
Then provided Vf ^ 0, 

P{f) = ^) y 

r Vf\r) 

Substituting into the rearranged Euler equation gives 



rvf 



ki 



rVf 



[17) 



(18) 



(19) 



where c is a function of p, and thus a function of Vr. This now completely specifies the 
force profile in terms of the desired velocity profile, Vf, v^, the equation of state, and a 
single integration constant ki. Calculating the fluid acceleration leads to 



V ■ V)v = < Vf{r)drVf{r) — 



which can be rearranged to yield 

1 



r + I Vf{r)drVQ{r) + 



Vr{r) v^{r) 



0, 



{v- V)v 



dr[Vr{ry] - 



21 f^eWM ^ , fMr) 



r + 



dr[r Vg{r)] \ e 



(20) 



(21) 



2 ' ■ ' r 

Finally, decomposing the external force into radial and tangential [torque-producing 
directions] we have 



ff = f-r^ki 



1 

rVf 



+ dr 



1 

rVf 



and 



The radial equation can undergo one further simplification to yield 
ff^kil- 



fi 


r c^i 


1 r 


1 - 

2 







drVf 



(22) 



(23) 



(24) 



Summarizing: 

• In the geometric acoustics regime the acoustic line-element for the most general 
[time-independent cylindrically symmetric collapsing/ expanding] line vortex is 

,2 / 2 ,.2 



ds^ oc - {c^ -Vf- vj) dt^ 



^ 2vf dr dt - 2vs r d9 dt 

2 , „2i/i2 , 1 2 



+ dr' + r'de'' + dz\ (25) 

We again reiterate that — given a barotropic equation of state p{p) — once the 
velocity profile Vf and Vg is specified, then up to a single integration constant ki, 
the density and speed of sound are no longer free but are fixed by the continuity 
equation and the equation of state respectively. Furthermore, the Euler equation 
then tells you exactly how much external force is required to set up the fluid flow. 
There are still two freely specifiable functions which we can take to be the two 
components of velocity. 
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• In the physical acoustics regime the acoustic hne-element is 




(26) 



The major difference for physical acoustics is that for technical reasons the massless 
curved-space Klein-Gordon equation [d'Alembertian wave equation] can only be 
derived if the flow has zero vorticity. This requires Wg(r) = k2/r and hence = 0, 
so that the flow would be un-torqued. More precisely, the d'Alembertian wave 
equation is a good approximation as long as the frequency of the wave is high 
compared to the vorticity. However, in the presence of significant torque and 
vorticity, a more complicated wave equation holds [7], but that wave equation 
requires additional geometrical structure beyond the effective metric, and so is 
not suitable for developing general relativity analogue models. (Though this more 
comphcated set of coupled PDEs is of direct physical interest in its own right.) § 

There are several special cases of particular interest: 

- No vorticity, V x v = 0; 

- No angular torque, /g = 0; 

- No radial flow v.f = Q\ 

- No radial force, /,-. = 0; 

- No external force, / = 0; 

which we will now explore in more detail. 

Zero vorticity/ zero torque: If we assume zero vorticity, Vxv = 0, the above calculation 
simplifies considerably, since then 

= 7 . (27) 

which implies f^ — 0. Conversely, if we assume zero torque, then (assuming ki ^ 0) the 
vorticity is zero. 

But note that the simple relationship zero torque <(=^ zero vorticity requires the 
assumption of nonzero radial flow. With zero radial flow the torque is always zero for a 
time independent flow, regardless of whether or not the flow is vorticity free. 

Zero radial force: Assuming zero radial force, ff = 0, and assuming ki ^ 0, one finds 



§ We shall subsequently see that the "equivalent Kerr vortex" is not irrotational — but the vorticity 
is proportional to the angular momentum, so there is a large parameter regime in which the effect of 
vorticity is negligible. 




(28) 
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Thus the radial and angular parts of the background velocity are now dependent on each 
other. Once you have chosen, e.g., '^e^'^) and c(r), a differential equation constrains vf. 

.21 c^ + vl 

^ . (29) 



1-^ 



Zero external force: If wc assume zero external force, / = 0, both the radial and 
angular external forces are zero. Then assuming fci 7^ 



1 - 



- + %■ (30) 



Since c depends on p via the barotropic equation, and p depends on Vr via the continuity 
equation, this is actually a rather complicated nonlinear ODE for Vr{r). (For zero radial 
flow this reduces to the tautology = and we must adopt the analysis of the previous 
subsection.) 

In general the ergosurface is defined by the location where the flow goes supersonic 
vl + v] = (31) 

while the horizon is defined by the location where 

vj = (32) 

Note that horizons can form in three rather different ways: 

• Vf = —c ^0 — a black hole horizon. 

• Vf. = +c 7^ — a white hole horizon. 

• Vf = ±c = — a bifurcate horizon. 

In analogue models it is most usual to keep c 7^ and use fluid flow to generate 
the horizon, this is the case for instance in the Painleve-Gullstrand version of the 
Schwarzschild geometry [4, 5]. The alternate possibility of letting both Vf and 
c ^ to obtain a horizon is not the most obvious construction from the point of 
view of acoustic geometries, but cannot a priori be excluded on either mathematical 
or physical grounds. Indeed, it is this less obvious manner of implementing the 
acoustic geometry that most closely parallels the analysis of Schwarzschild black holes 
in curvature coordinates or isotropic coordinates, and we shall soon see that this route 
is preferred when investigating the Kerr spacetime. 

3. The role of dimension 

The role of spacetime dimension in these acoustic geometries is sometimes a bit 
surprising and potentially confusing. This is important because there is a real physical 
distinction between truly two-dimensional systems and effectively two-dimensional 
systems in the form of three-dimensional systems with cylindrical symmetry. We 
emphasise that in cartesian coordinates the wave equation 

^ ^rilz A=o, (33) 



dxi^ \ dx^ 
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where 



-p v'^/c^ I p {S^^ — v'^v^/c^} 



(34) 



holds independent of the dimensionahty of spacetime. It depends only on the Euler 
equation, the continuity equation, a barotropic equation of state, and the assumption 
of irrotational flow [3] . 

Introducing the acoustic metric g'^'^, defined by 

1 



-9 9 



1X1/ . 



9 = 



(35) 



dct(^^'^) 

the wave equation (33) corresponds to the massless Klein-Gordon equation 
[d'Alembertian wave equation] in a curved space-time with contravariant metric tensor: 



2/(d-l) 


" -1/C2 






— v/c^ 





(36) 



where d is the dimension of space (not spacetime). 
The covariant acoustic metric is 



9fii/ 



& 



1) 


" - (c^ - v^) 


->T 
—V 




— * 

—V 





(37) 



3.1. d = 3 

The acoustic line-element for three space and one time dimension reads 



9n2/ 



- (c2 - V^) 


-'T 
—V 


— * 

—V 


IsxS 



(38) 



This is the primary case of interest in this article. 



3.2. d = 2 

The acoustic line-element for two space and one time dimension reads 



- (C2 - V^) 


1 -'T 
—V 


—V 


1 I2x2 



(39) 



This situation would be appropriate when dealing with surface waves or excitations 
confined to a particular substrate, for example at the surface between two superfluids. 
An important physical point is that, due to the fact that one can always find a conformal 
transformation to make the two-dimensional spatial slice flat [20] , essentially all possible 
2+1 metrics can in principle be reproduced by acoustic metrics [21]. (The only real 
restriction is the quite physical constraint that there be no closed timelike curves in the 
spacetime.) 
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3.3. d^l 

The naive form of the acoustic metric in (1+1) dimensions is ill-defined, because the 
conformal factor is raised to a formally infinite power — this is a side effect of the 
well-known conformal invariance of the Laplacian in 2 dimensions. The wave equation 
in terms of f^^ continues to make good sense — it is only the step from f^^^ to the 
effective metric that breaks down. 

Note that this issue only presents a difficulty for physical systems that are 
intrinsically one-dimensional. A three-dimensional system with plane symmetry, or 
a two-dimensional system with line symmetry, provides a perfectly well behaved model 
for (1-1-1) dimensions, as in the cases d — ?> and d — 2 above. 



4. The Kerr equator 



To compare the vortex acoustic geometry to the physical Kerr geometry of a rotating 
black hole [22], consider the equatorial slice 6* = 7r/2 in Boyer-Lindquist coordinates [23]: 

2m , , , , , o dr^ 



(ds^)(2+i) = -dt' + — (dt - a d(t>Y + / ^ 2 / 2 

r 1 — 2m/r + a'^ /r'^ 

We would like to put this into the form of an "acoustic metric" 

_{c2 _ h^n 



+ (r^ + a^) d0^ 



h. 



If we look at the 2-d r-0 plane, the metric is 

dr^ 



ids') 



(2) 



+ + a" + 



2ma^ 



d<f)'. 



(40) 



(41) 



(42) 



1 — 2m/r + a? /r"^ 

Now it is well-known that any 2-d geometry is locally conformally flat [20] , though this 
fact is certainly not manifest in these particular coordinates. Introduce a new radial 
coordinate f such that: 



dr^ 



2m/r + /r"^ 



This implies 



and 



+ + 



dr^ 



2ma^ 



1 — 2m/r + a'^/r'^ 
leading to the differential equation 
1 df(r) 



r'' + a'^ + 



= Q(f)' df^ 



2ma^ 



d0' = Q(f)2 [df2 + Pd02]. (43) 



(44) 



r{r) dr - 2m/ r + a^/r'^^r'^ + + 2ma^/r 



(45) 



(46) 
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— 2m/r + jr'^^Jr'^ + + 2ma^ /i 
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(47) 



The normalization is most easily fixed by considering the m = = a case, in which case 
f = r, and then using this to write the general case as 



f(r) = r exp 



1 



1 



dr 



(48) 



— 2m/r + /r'^^Jr'^ + a? + 2ma^ jr r 

where f is simply a dummy variable of integration. If a = this integral can be 
performed in terms of elementary functions 

dr 



'a^I — 2m/r 



~ In — 2mr + r — [r > 2m], 



so that 



r r 



- (^V — 2mr + r — , 



(49) 



(50) 



though for m ^ and general a no simple analytic form holds. Similarly, for m — and 

a 7^ it is easy to show that 

f(r) = Vr^ + a^, (51) 

though for a 7^ and general m no simple analytic form holds. 

Nevertheless, since we have an exact [if formal] expression for f(r) we can formally 
invert it to assert the existence of a function r(f). It is most useful to write f = r F{r), 
with limr^oo -^('") = 1, and to write the inverse function as r = f H{r) with the 
corresponding limit lim^^^oo -f^('^) = 1- Even if we cannot write explicit closed form 
expressions for F{r) and H{r) there is no difficulty in calculating them numerically, 
or in developing series expansions for these quantities, or even in developing graphical 
representations. 

We now evaluate the conformal factor as 



+ + 2ma?/r 



H r 



1 + 



2ma^ 



with r considered as a function of f, which now yields 



2m 



(ds2)(2+i) = -dt^ + — {dt^ - 2a d(f) dt) + n{ry [dP + P d^^]. 



(52) 



(53) 



Equivalently 



(d.^)(2+i) = n(f)^ -0(f) 



^ 2m 
r 



dt' 



^{r) d(j) dt 



+ [df^ + d(p 



(54) 



This now lets us pick off the coefficients of the equivalent acoustic metric. For the overall 
conformal factor 



^ = Q(f)^ = //^(f) 1 + ^ + 



a 



(55) 
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For the azimuthal "flow" 



^;^-2 



2am 2am 



r r 

In terms of orthonormal components 

2am 



H-\f) 1 + + 



2ma^ 



E-\f) 1 + 3 + 



2ma' 



(56) 



(57) 



This is, as expected, a vortex geometry. FinaUy for the "coordinate speed of hght", 
corresponding to the speed of sound in the analogue geometry 



= Q(f) 



-2 



2m 
r 



^ 2m 




r 


+ ,~2 ,2 1 



Q(f)-* <^ Q(f) 



The speed of sound can be rearranged a httle 

2ma^ 



1 + 



1 - 



2m 



4a^m^ 



+ 



This can now be further simphfied to obtain 
and finally leads to 



TT/~\ /-, 2m 
p^H{r) \1 + -, 



(58) 
(59) 

(60) 

(61) 

(62) 



with r implicitly a function of r. Note that the velociy field of the "equivalent Kerr 
vortex" is not irrotational. Application of Stokes' theorem quickly yields 



1 d 



27rr dr 



1 d 



1 d 



r dr 



r df 



2amr 



7.2 



1 



a 2ma" 



2\ -1 



which we can somewhat simplify to 



4am 



2ma^ 
1 + ^ + 



2\ -1 



(63) 



(64) 



+ fl r a/I — 2m/r + a^/r 



9_ 
9r 



2 am 



a 2ma" 
1 + 7? + ^ 



Note that (as should be expected), the vorticity is proportional to an overall factor of a, 
and so is proportional to the angular momentum of the Kerr geometry we are ultimately 
interested in. 

While f is the radial coordinate in which the space part of the acoustic geometry 
is conformally fiat, so that f is the "physical" radial coordinate that corresponds to 
distances measured in the laboratory where the vortex has been set up, this particular 
radial coordinate is also mathematically rather difficult to work with. For some purposes 
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it is more useful to present the coefficients of the acoustic metric as functions of r, using 
the relationship F{r) — 1/H{f). Then we have: 



p(r) 
c(r) 



2am 



a 2ma^ 

1 + ;:^ + ^ 



F(r) 



a 2ma^ 

1 + ;:^ + ^ 



-1 



2m ^, , 
+ — Fir) 



c{r) ^ \/l 

Finally this yields 

p(r) = F-(r)yi ^ 
and the explicit if slightly complicated result that 



a 2ma^ 
1 + ^ + ^ 



2m 



r.2 ' 



1 



- S dr 



(65) 
(66) 
(67) 

(68) 
(69) 



F{r) — exp . , _ _ _ 

y 1 — 2m/r + a^/r^ y/r^ + + 2ma^/r r 

One of the advantages of writing things this way, as functions of r, is that it is now 
simple to find the locations of the horizon and ergosphere. 

The ergo-surface is defined by — — 0, equivalent to the vanishing of the gu 
component of the metric. This occurs at 

rE = 2m. (70) 

The horizon is determined by the vanishing of c, (recall that Vf = 0), this requires 
solving a simple quadratic with the result 

rn = Tn + \/m^ — < rE- (71) 

These results agree, as of course they must, with standard known results for the Kerr 
metric. It is easy to check that F{rH) is finite (the function inside the exponential is 
integrable as long as the Kerr geometry is non-extremal), and that F{rE) is finite. This 
then provides some simple consistency checks on the geometry: 

Pirn) = = cim); (72) 



F{rH) 



c[rH); 
finite; 

G{rH)~^ = finite; 
a 



(73) 
(74) 
(75) 



cu(rfl-) 
c{rE) 



(76) 
(77) 



2mrH 
= \-^4>{rE)\- 

In particular if we compare the vorticity at the horizon, uj{rH) — (i/{2mrH), with the 
surface gravity 

2m 2a^ 
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we find a large range of parameters for which the peak frequency of the Hawking 
spectrum is high compared to the frequency scale set by the vorticity — this acts 
to suppress any complications specifically coming from the vorticity. (Of course in the 
extremal limit, where Th ^ the vorticity at the horizon becomes infinite. In the 
extremal limit one would need to make a more careful analysis of the effect of nonzero 
vorticity on the Hawking radiation.) 

Prom the point of view of the acoustic analogue the region inside the horizon is 
unphysical. In the unphysical region the density of the fluid is zero and the concept of 
sound meaningless. In the physical region outside the horizon the flow has zero radial 
velocity, and zero torque, but is not irrotational. Since by fitting the equatorial slice of 
Kerr to a generic acoustic geometry we have fixed p(r) and (?{r) as functions of r [and 
so also as functions of f] it follows that p(r) is no longer free, but is instead determined 
by the geometry. From there, we see that the EOS p{p) is determined, as is the external 
force ff{r). The net result is that we can [in principle] simulate the Kerr equator exactly, 
but at the cost of a very specific fine-tuning of both the equation of state and the 
external force fr{r). 

In figure 3 we present an illustrative graph of the speed of sound and velocity of 
the fiuid. Key points to notice are that the speed of sound (brown curve) and velocity 
of fiow (black curve) intersect on the ergosurface (green line), that the speed of sound 
goes to zero at the horizon (blue line), and that the velocity of fluid flow is flnite on 
the horizon. Finally the inner horizon is represented by the red vertical line. While 
the region inside the outer horizon is not physically meaningful, one can nevertheless 
mathematically extend some (not all) of the features of the flow inside the horizon. 
This is another example of the fact that analytic continuation across horizons, which 
is a standard tool in general relativity, can sometimes fail for physical reasons when 
one considers Lorentzian geometries based on physical models that differ from standard 
general relativity. This point is more fully discussed in [6, 24, 25]. 

In flgiirc 4 wc present an illustrative graph of the fluid density. Note that the density 
asymptotes to a constant at large radius, and approaches zero at the horizon. Finally, in 
flgure 5 we present an illustrative graph of the conformal factor Vt^ which remains flnite 
at the horizon and ergosurface, and asymptotes to unity at large distances from the 
vortex core. Figures 3-5 all correspond to m = 2 and a = 1, and have been calculated 
using a 30'th-order Taylor series expansion. || 

5. Discussion 

We have shown that the Kerr equator can [in principle] be exactly simulated by an 
acoustic analogue based on a vortex flow with a very speciflc equation of state and 
subjected to a very speciflc external force. Furthermore we have as a result of the analysis 
also seen that such an analogue would have to be very speciflcally and deliberately 
engineered. Thus the results of this investigation are to some extent mixed, and are 

II A suitable Maple worksheet is available from the authors. 
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Figure 3. Illustrative graph of the speed of sound and rotational velocity of the vortex 
as a function of r. (m — 2, a — 1, 30'th-order Taylor series.) 
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Figure 4. Illustrative graph of the Quid density as a function of r. (m = 2, a ^ 1, 
30'th-order Taylor series.) 

more useful for theoretical investigations, and for the gaining of insight into the nature 
of the Kerr geometry, than they are for actual laboratory construction of a vortex 
simulating the Kerr equator. 

One of the technical surprises of the analysis was that the Doran [26, 27] coordinates 
(the natural generalization of the Painleve-Gullstrand coordinates that worked so well 
for the Schwarzschild geometry) did not lead to a useful acoustic metric even on the 
equator of the Kerr spacetime. Ultimately, this can be traced back to the fact that 
in Doran coordinates the space part of the Kerr metric is non-diagonal, even on the 
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Figure 5. Illustrative graph of the conformal factor U as a function of r. (m — 2, 
a = 1, 30'th-order Taylor series.) 



equator, and that no simple coordinate change can remove the off-diagonal elements. In 
Boyer-Lindquist coordinates however the space part of the metric is at least diagonal, 
and the f coordinate introduced above makes the spatial part of the equatorial geometry 
conformally flat. This f coordinate is thus closely related to the radial coordinate in the 
isotropic version of the Schwarzschild solution and f in fact reduces to that isotropic 
radial coordinate as a — 0. 

However, in a more general physical setting the Doran [26, 27] coordinates might 
still be useful for extending the acoustic analogue away from the equator, and into the 
bulk of the domain of outer communication. To do this we would have to extend and 
modify the notion of acoustic geometry. The fact that the spatial slices of the Kerr 
geometry are never conformally fiat [17, 18, 19] forces any attempt at extending the 
acoustic analogy to consider a more general class of acoustic spacetimes. A more general 
physical context that may prove suitable in this regard is anisotropic fluid media, and we 
are currently investigating this possibility. The simplest anisotropic fluid media are the 
classical liquid crystals [28, 29], which in the present context suffer from the defect that 
they possess considerable viscosity and exist chiefly at room temperature — this makes 
then unsuitable for the construction of analogue horizons, and particularly unsuitable 
for the investigation of quantum aspects of analogue horizons. Much more promising in 
this regard are the quantum liquid crystals. These are anisotropic superfluids such as, in 
particular, 3He-A [21, 30]. The anisotropic superfluids are generally characterized by the 
presence of a non-scalar order parameter; they behave as superfluid liquid crystals, with 
zero friction at T = 0. The premier example of an anisotropic superfluid is fermionic 
3He-A, where the effective gravity for fermions is described by vierbeins, though several 
types of Bose and Fermi condensates of cold atoms also exhibit similar behaviour [21, 30]. 

In the system we have presented above we have found that two of the main features 
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of the Kerr geometry appear: the horizon (outer horizon) and ergo-region. By shght 
modification of standard arguments, these geometric features are expected to lead to 
the quantum phenomena of super-radiance and Hawking radiation [1, 2, 3, 4, 5]. The 
super-radiance corresponds to the reflection and amphflcation of a wave packet at the 
ergo surface [31]. (See also [32, 33].) The model, because it contains a horizon, also 
satisfies the basic requirements for the existence of Hawking radiation, which would 
now be a thermal bath of phonons emitted from the horizon [1]. ^ These are two of 
the most fundamental physics reasons for being interested in analogue models [5]. A 
subtlety in the argument that leads to Hawking radiation arises from the way that the 
choice of coordinates seems to infiuence the choice of quantum vacuum state. In the 
presence of a horizon, the choice of quantum vacuum is no longer unique, and standard 
choices for the quantum vacuum are the Unruh, Hartle-Hawking, and Boulware states. 
It is the Unruh vacuum that corresponds (for either static or stationary black holes) 
to Hawking radiation into an otherwise empty spacetime, while the Hartle-Hawking 
vacuum (for a static black hole) corresponds to a black hole in thermal equilibrium with 
its environment (at the Hawking temperature). For a stationary [non-static] black hole 
the Hartle-Hawking vacuum state does not strictly speaking exist [34], but there are 
quasi-Hartle-Hawking quantum states that possess most of the relevant features [35] . In 
general relativity, because physics is coordinate independent, the choice of vacuum state 
is manifestly independent of the choice of coordinate system. In the analogue spacetimes 
considered in this article, because the preferred choice of coordinate system is intimately 
tied to the flowing medium, the situation is perhaps less clear. The coordinate system 
we have adopted is invariant under combined time reversal and parity, which might 
tempt one to feel that a pseudo-Hartle-Hawking vacuum is the most natural one. On 
the other hand, if we think of building up the vortex from a fluid that is initially at 
rest, then there is clearly a preferred time direction, and the "white hole" component 
of the maximally extended horizon would exist only as a mathematical artifact. In this 
physical situation the Unruh vacuum is the most natural one. We feel that there is an 
issue worth investigating here, to which we hope to return in some future article. 

Finally, since astrophysically all black holes are expected to exhibit some degree 
of rotation, it is clear that an understanding of the influence of rotation on analogue 
models is important if one wishes to connect the analogue gravity programme back 
to astrophysical observations. In conclusion, there are a number of basic physics 
reasons for being interested in acoustic analogues of the Kerr geometry, and a number 
of interesting directions in which the present analysis might be extended. 

% Furthermore, there is a strong feeling that the cosmological constant problem in elementary particle 

physics might be related to a mis-identification of the fundamental degrees of freedom [21, 30]. 

+ An interesting attempt in the opposite direction [36, 37] is the recent work that analyzes accretion 

flow onto a black hole in terms of a dumb hole superimposed upon a general relativity black hole. 
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Appendix: 

Isotropic version of the Schwarzschild geometry cast in acoustic form 

The isotropic version of the Schwarzschild geometry can be cast exactly into acoustic 
form, without any extraneous conformal factors. Since this is not the standard way of 
viewing the acoustic analogue of the Schwarzschild geometry [4], and since the algebra 
is simple enough to be done in closed form, it is worthwhile taking a small detour. 
In isotropic coordinates the Schwarzschild geometry reads 

= + + ?^)' i^^' + + ^i^' ^ ^'^')} ' (^9) 

and so in these coordinates the acoustic analogue corresponds to 

V =0; (80) 

P=Poo(i-^); (81) 



The external force required to hold this configuration in place against the pressure 
gradient is 

" dr " dr " (1 + f )6 2r3- 

The pressure itself (normalizing poo 0) is then 

P{r) = -p^c^m ^ 480 (1 + f )^ j' ^''^ 

and by eliminating r in favour of p using 

m 



(85) 



~ pIPoo 

we can deduce the equation of state 
/X / .2 /^20 - hp/ Poo - 4^1 - p/poo - p/poo\/l- pIPoo , i^p. 

At the horizon, which occurs at = m/2 in these coordinates, both p — and c — 
in the acoustic analogue, while 

PH = (87) 

is finite, so there is a finite pressure drop between asymptotic infinity and the horizon. 
Everything is now simple enough to be fully explicit, and although we now have an 
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analogue model that reproduces the exterior part of the Schwarzschild geometry exactly^ 
we can also clearly see the two forms in which fine tuning arises — in specifying the 
external force, and in the equation of state. 

More generally, it is a standard result that any spherically symmetric geometry can 
be put into isotropic form 

ds^ = - exp[-2*] dt^ + exp[-2*] {dr^ + r'^{de'^ + sin^ 9 dc/)^)} . (88) 

This can be put into acoustic form in a particularly simple manner by setting 

v{r) = 0; (89) 
p(r) =poo exp[-$-*]; (90) 
c(r) =Coo expf^'-^]. (91) 

The pressure can now be formally evaluated as 

/oo 
exp[*-3$]9,($ + *)dr, (92) 

and comparison with p{r) can be used to construct a formal equation of state, p{p). 
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